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O |' We discuss the theoretical treatment of non-global observables, those quantities that are sensi- 

tive only to radiation in a restricted region of phase space, and describe how large 'non-global' 
logarithms arise when we veto the energy flowing into the restricted region. The phenomeno- 
| logical impact of non-global logarithms is then discussed, drawing on examples from event 

shapes in DIS and energy-flow observables in 2-jet systems. We then describe techniques to 
reduce the numerical importance of non-global logarithms, looking at clustering algorithms in 
energy flow observables and the study of associated distribution of multiple observables. 
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1 Introduction and Theory 

Recently a distinction has been introduced between so-called global and non-global QCD ob- 
servables.^ The former are sensitive to emissions in all directions, while the latter are sensitive 
only to emissions in some restricted angular region, for an example a jet or a hemisphere. Ob- 
vious examples of non-global (NG) observables are properties of individual jets (invariant mass, 
numbers of subjets). Many other common QCD observables are also non-global, including def- 
initions of diffraction based on rapidity gaps (whether in terms of particles or energy flow); 
isolation criteria for photons (or other particles); distributions of interjet energy flow; or even 
the original Sterman- Weinberg jet definitional 

The question of globalness becomes of particular relevance whenever one places a severe 
restriction on the energy E (or in some cases, transverse energy) flowing into the observed 
region. In such a situation the perturbative series develops logarithmically enhanced terms at 
all orders, at the very least single logs a™ ln n E/Q, where Q is the hard scale. For E <C Q such 
a series needs to be resummed. For global observables it has been showrPUhat the resummation 
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Figure 1: Real and virtual contributions to 
the emission of gluon (2) from a qqg system. 

in other words one can make an 



can be carried out to single logarithmic (SL) accuracy, essentially by using the approximation 
of independent emission off the underlying Born event. 

Until recently it had universally been assumecp'that 
this approximation could be used more generally. How- 
ever it turns out that for non-global observables, SL re- 
summation is more subtle. ^ This is illustrated in fig.Q 
which shows left and right hemispheres (He, Hn) of a 
2-jet e + e~ event, and two emissions going into oppo- 
site hemispheres such that Q » Bi > i?2. A global 
observable would for example measure the sum of the 
two gluon energies. Since E\ 3> Ei, placing a restric- 
tion E max on the sum is equivalent to placing it directly 
on Ei and there is cancellation between the real pro- 
duction and virtual loop contribution for gluon 2. Be- 
cause of this cancellation, one is free to 'mistreat' the 
way gluon 2 is emitted and pretend it is emitted from 
the simpler qq system, ignoring the presence of gluon 1 
independent emission approximation. 

Now suppose we have an observable that measures the energy only in H-jz- The limit is placed 
just on gluon 2, Ei < E max . On the other hand the loop contribution has an effective limit 
-E^virtuai ^ E\ and the mismatch between these two limits leads to a logarithmic enhancement 
lnE±/E2- After integrating over E\ one finds an overall contribution a 2 s In 2 Q/E max . Making an 
independent emission approximation, one would obtain the wrong coefficient for this term. The 
difference between the true answer (based on the coherent emission of gluon 2 from the qqgi 
system) and the independent emission result is termed a 'non-global logarithm' (NGL). 

At this two-gluon level, non-global logarithms are 
essentially an edge effect: it is only close to the bound- 
ary between measurement and non-measurement that 
one is sensitive to the difference between independent 
emission and the true two-gluon emission pattern. This 
is illustrated in figure[2]which shows (through the colour 
shading) the contribution to the NGL as a function of 
the two gluons' rapidities rj and azimuths 4> for the case 
in which the 'measurement' is carried out in a square 
patch |?7| < 1, |0| < l. 6 One sees clearly that the largest 
contribution to the NG term is concentrated around the 
borders of the patch. 

In the region where a s lnQ/E max is of order 1, one 
needs to understand non-global effects at all orders 
a™ In™ Q/E max . This amounts to accounting for the 
coherent radiation into the observed region of soft glu- 
ons from arbitrarily complex ensembles of harder (but 
still energy-ordered) gluons in the non-observed region. 
It turns out that given a naive resummed calculation 
based on an independent emission picture, non-global 
effects can be accounted for by a multiplicative correction factor S(t), where t is the running- 
coupling generalisation of & InQ / E max : 

dE t a s (E t ^ 




Figure 2: Contribution to the second or- 
der NGL for measurement in a square patch 
shown as a function of rapidity and azimuth 
Darker shades correspond to larger contribu- 
tions. 
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When viewed in colour, the blue and red shadings are therefore respectively for the unmeasured and measured 
gluons. 
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Figure 3: The function S(t) for different geometries of 
observed region in a two-jet event. Slices are defined 
as containing all azimuths but a restricted range of 
rapidity. 





Figure 4: Above: forbidding coherent radiation from 
an unobserved region (grey gluons) into the observed 
region (black gluons with crosses); below: schematic 
representation of the actual distribution of gluons in r\ 
and <f> (observed region is hashed). 



It is sometimes useful also to write the expansion of S, S(t) = ^ l ■> an d f° r example the 

two-gluon contribution discussed above gives us £2- The higher order terms are more difficult 
to calculate because of the complicated colour structures that appear for emission from multi- 
gluon configurations, and also because of the geometry. The colour problem has so far only been 
partially solved: using the large-iV c approximation one can restrict one's attention to planar 
graphs, equivalent to considering emission from sets of independent colour dipoles0 

Two equivalent approaches have been proposed to deal with the geometry dependence. In 
practice, the simplest way of solving the problem of the geometry dependence is through a Monte 
Carlo branching algorithm, ^ in which the original qq dipole branches into two dipoles qg and 
gq. As one increases the logarithm t each new dipole can itself branch and one iteratively builds 
up an ensemble of energy-ordered gluons with the correct (large- N c ) angular distribution. This 
algorithm is similar to that of the Ariadne event generator!^ One then determines S(t) by taking 
the number of events at scale t that are free of emissions in the observed region and dividing it 
by the number of events that would have been expected on the basis of an independent emission 
picture. 

Results for S{t) are shown in fig. for various geometries of observed regions. From a theo- 
retical point of view the most remarkable feature of these curves is that modulo normalisation, 
they all have the same t-dependence (in contrast, the fixed order £2 terms differ by more than a 
factor of two between the different geometries). The explanation proposed for this observation^ 
was so the so-called 'buffer mechanism', illustrated in fig. ^ the hypothesis is that the easiest 
way of forbidding 'secondary' emissions into the observed region is actually to forbid primary 
emissions close to the observed region. Accordingly for intermediate scales if < t one finds a 
buffer region, free of emissions, surrounding the observed region. Larger differences t — t' imply a 
larger buffer region. As a result the dynamics governing the large-t behaviour of S mostly occurs 
far from the observed the region and is not affected by the geometry of the observed region. 
Assuming that the buffer region is of size 5r/b u ffcr — cC^i — t') one comes to the conclusion that 



5(i)~ e - 2cC ^ 2 . 



(2) 



These arguments were placed on a mathematically sound footing by Banfi, Marchesini and 
SmyeP. Firstly they introduced an alternative but equivalent approach to the problem of the 
geometry, in terms of a non-linear integral equation: 



dt S ab (t) = [ ^ w ab {k) (uf b \{t) S ak (t) S kb (t) - S ab (t)) , (3) 

J unmeasured 

where S ab {t) is the non-global correction factor for a dipole ab, w ab (k) is the weight for emission 
of a soft gluon k from ab and U^ k (t) accounts for the different 'primary' emission contributions 
for dipoles ak, kb compared to dipole ab. They were then able to demonstrate the existence 
of scaling solutions to (J3J), proving the buffer mechanism, with c evaluated numerically to be 
2.5 ±0.25%. They also evaluated a number of subleading corrections. They pointed out however 
that in practice (J2J) applies only for t > 0.5, i.e. very asymptotic values oft. Phenomenologically, 
t is limited to be < 0.2 and one should use the full solutions to S. 

1.1 Warnings for practitioners, a.k.a. Zoology 

The definition of a non-global observable, given above, is one that is sensitive only to emissions 
(from the Born event) in a restricted angular region. However the situations in which NGLs can 
appear are rather subtle. 

Firstly there exist observables that measure only a subset of the particles, but which never- 
theless are global. An example is the heavy-hemisphere invariant squared mass in e + e~ . Though 
only one hemisphere is measured, the observable is global because it is always the heavier hemi- 
sphere that is measured — in a situation such as fig. ^ the heavier hemisphere is always the one 
with the harder gluon (1) and one therefore never has the situation of a harder unobserved par- 
ticle radiating into the observed region. Other examples include certain DIS event shapes (e.g. 
B z e, t z q) where the measurement is carried out in one current hemisphere, but for which con- 
servation of momentum causes an indirect sensitivity to emissions in the unobserved hemisphere. 
Such observables are known as indirectly global. 

Another subtle case is that of observables referred to as dis continuously global. Such observ- 
ables have different parametric sensitivities to emissions in different directions, e.g. v ~ Ef/Q 2 
for emissions in one hemisphere, v ~ Et/Q for the other. Placing a limit on v corresponds to 
different limits on E% in the two hemisphere (e.g. \fvQ and v Q respectively) and one finds NGLs 
as for simple non-global observables, but with the appropriate replacement of the integration 
limits for t in eq. (|T|)!^^ 

The most subtle situation is perhaps that of dynamically discontinuously global observables. 
These are observables which for a single emission appear global (typically indirectly global). 
However they involve non-linearities such that for the configurations of emissions that are most 
common (e.g. given a certain value of the observable) they develop different effective paramet- 
ric dependence on emissions in different regions.^ One example of such an observable is the 
broadening B z e in DIS. 

2 Phenomenological implications 

In this section we will examine the phenomenological impact of NGLs in a variety of non-global 
observables. We will look at the invariant-squared jet mass in (although the conclusions 

we draw will apply to e + e~ non-global event shapes as well) and at 2-jet energy flows.^^ In 
both kinds of observable we will see that by neglecting the non-global logarithm suppression 
factor, one overestimates the distributions and cross sections by a considerable, and certainly 
phenomenologically relevant, amount. 



2.1 NGLs in DIS and e + e event shapes 

Event shapes (ES), and their distributions,^ are widely used to compare the predictions of 
perturbative QCD with experimental observation. As an example of a non-global event shape, 
we will use the DIS invariant squared jet mass p, defined as 

4(s Wc |Pi|) 

where the summation is over all particles in the current-hemisphere (Ti-c)- When calculating 
the cross section for p to be smaller than some given value, one finds that in the exclusive limit, 
small p, each power of the coupling is accompanied by up to two powers of In p, associated with 
with soft and collinear divergences. These terms need to be resummed. The standard accuracy 
is next-to-leading logarithmic, which in this case implies the inclusion of all single-logarithmic 
terms, i.e. the same accuracy as NG contributions. 

In figure El we show the effect of the non-global logarithms on the resummed distribution 
for p. The broken line is the resummation without NGLs, and the solid line shows the result once 
they are included through the function S discussed above. Neglecting the non-global logarithms 
leads to an overestimation of the peak- height by around 40%, while the effect is smaller in other 
parts of the distribution. In practice, the resummed calculation is usually matched to a fixed- 
order calculation, which includes the NGLs to order a 2 s and so reduces the impact of neglecting 
the NGLs in the final result. 



(4) 



2.2 NGLs in 2 jet energy flow observables 

Let us now examine the impact of non-global logarithms on energy flow observables.^ Such 
observables have attracted con siderable interest in recent time as an infrared-safe wa y of studying 
gaps-between-jets processe^^ and the underlying event in hadron-hadron collisions!^ We shall 
take as our observable the total amount of transverse energy flowing into a restricted region of 
phase space fi. The probability for the amount of transverse energy to be smaller than some 
value Qn is 

£n(Qn) = -/ dE t —. (5) 
As discussed in section ^ this can be written as the product of two contributions, 



En(Qn)=S(t)xE P (t), (6) 

where the function Ep is based on an independent gluon emission approximation ('primary' 
radiation into while S(t) accounts for secondary coherently radiated emission into f2 and t is 
defined as the integral of a s /2ir between Qq and Q. We can understand the phenomenological 
implications of the non-global logarithms on this observable by comparing the result for En, 
calculated with only primary emissions, with the full result for En, which accounts also for non- 
global logarithms. We do this in figure El where f2 is a slice in rapidity of width At] = 1.0. The 
plot clearly show the phenomenological impact of the NGLs on this observable; the suppression 
is significant, particularly at larger values of the effective logarithm t. The typical energies of 
current colliders corresponds to about t = 0.15 and if we take this as our reference value, the 
inclusion of non-global logarithms increases the suppression, relative to the primary-only case, 
by a factor of about 1.65. Similar results are found for other definitions of f2, for example a 
patch in phase space bounded in rapidity and azimuthal angle. 
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Figure 5: The impact of non-global logarithms (which 
enter at the NLL level) on a non-global DIS event 
shape. 
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Figure 6: The impact of non-global logarithms on the 
2-jet energy flow observable Eq. The upper plot shows 
the relation between Qn, Q and t. 



3 Controlling non-global logarithms 

In this section we will describe ways of controlling, or taming, non-global logarithms. Such a 
study is important given that our methods for resumming NG logarithms are both approximate 
(large- Nq limit) and cumbersome (numerical or asymptotic). We will look at two different 
approaches: minimising the numerical effect of the NGLs by clustering the final state, and also 
by examining associated distributions of pairs of observables. 

3.1 Clustering algorithms 

The applicatio n o f clustering algorithms'^ to the final state is motivated by recent HI and 
ZEUS analyses^! of gaps-between-jets processes at HERA. In these analyses, the inclusive k t 
algorithm is used to define the hadronic final state, and hence the rapidity gap Q, and a gap 
event is defined by a restriction of the total transverse energy into f2. This observable, sensitive 
to soft radiation into f2 only, is clearly non-global, and so sensitive to non-global logarithmic 
effects. However, as we will sho w, t he clustering procedure reduces the numerical importance 
of the NGLs for this observable,'^ and we can see why by looking at how the kt clustering 
algorithm works. The essential feature is that in an iterative algorithm over all the final state 
particles, it merges particles of lower transverse momentum into particles with higher transverse 
momentum, to produce pseudo-particle or mini-jets. Consider applying this algorithm to two 
gluons with strongly ordered transverse momentum, 

-Et,i ^ -Et,2> CO 

where gluon one, directed outside of the region $7, then in turn radiates gluon two into the region 
O. The strong ordering ensures this configuration produces a non-global logarithm. When we 
cluster this system, gluon two will be clustered into, or merged with, gluon one (and hence 
pulled out of the gap) if the two gluons are sufficiently close in the (rj, 4>) plane, 



(vi-V2) 2 + (0i-fo) 2 <R 2 , 



(8) 




Figure 7: The leading, S2 contribution to the non- 
global logarithm function S(t), with and without clus- 
tering. The saturation of S2 at high Aj; is seen for 
both the non-clustered and the clustered case, with a 
lower saturation value for the latter. 



Figure 8: The phenomenological impact of the kt clus- 
tering procedure is to reduce the numerical importance 
of non-global logarithms on this 2-jet energy flow ob- 
servable. 



where R is a parameter playing the role of a radius in the algorithm. Therefore to get a 
kinematical configuration which produces a NGL, the two gluons need to be sufficiently separated 
in the (77, <j>) plane to avoid being merged. Hence the kt algorithm 'cleans up' the restricted 
region of phase space, Q, and pulls soft gluons out of the gap. Figure [7| shows the leading, 
1S2 contribution (order af) to the NGL function S(t), without clustering (solid line) and with 
clustering (broken line.) In both cases S2 rapidly saturates at high Arj, which results from the 
fact that the NGLs are an edge effect, and the saturation value for the clustered case is smaller 
than that of the non-clustered case. This follows from the fact that, although the clustering 
algorithm pulls gluons out of the gap, gluons can still be sufficiently separated in the (77, cj>) 
plane to survive clustering and give a significant NGL contribution. Figure El shows the all- 
orders calculation of the full function E, where f2 is taken to be a slice in rapidity of A77 = 1.0. 
These figures allow us to see the phenomenological impact of the non-global logarithms on these 
observables, as they show the function Ep (only primary logarithms) and the full function E 
(primary and non-global logarithms) with and without clustering. We recall that t is around 
0.15 to 0.2 at the energies of current colliders, and as before we will take t = 0.15 as our reference 
value. The plot shows that the non-global logarithms cause a considerable suppression of E, 
relative to the primary-only result, and that by clustering the final state this suppression is 
reduced. At t = 0.15 the full result without clustering is suppressed relative to the primary only 
result by 1.65, and this suppression is reduced with clustering to around 1.2. 



3.2 Event shape/Energy flow Correlations 

Another way of reducing the phenomenological impact of NGLs that has been proposed is 
the study of associated distributions in two variables. In this work, one combines measurement 
of a jet shape V in the whole of phase space (for example thrust, V = 1 — T) and that of the 
transverse away-from-jets energy flow E OVL f The former is a global measurement and the latter 
is a non-global measurement. If the observable V selects 2-jet-like configurations, one measures 



the associated distribution 



^2ng(Q, V, E out ) 



(9) 



where Q is the hard scale. It has been shown that this distribution factorizes 



9 



s 2ng (Q, v, Eout) = e(q, v) ■ s out (yg, E out ) 



(10) 



where V) is the standard global distribution of V and T, out (VQ, E out ) contains the logarith- 
mic distribution in E out . This latter distribution, containing non-global logarithms is evaluated 
at the reduced scale VQ, and hen ce the logarithmic terms will be (a s log(V Q / E ont )) n • The work 
of Berger, Kucs and Stermar^^ considered the region in which VQ and E out were comparable, 
so that the NGLs give a negligible contribution. Thus, for a restricted subset of appropriately 
selected events, it is possible, to 'tune out' the non-global logarithmically enhanced terms in 
associated distributions. 

4 Conclusion 

To summarise, non-global logarithms are recently discovered contributions that arise in the 
distributions of any QCD observable sensitive only to emissions in a restricted part of phase 
space. They are phenomenologically important and significant progress has been made in re- 
summing them to all orders in the large- iV c limit. One of the main directions of current work 
focuses on understanding ways of designing observables so as to reduce the impact of non-global 
contributions. 
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